Four-body structure of A Li and AiV spin-dependent interaction 
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Two spin-doublet states of in A Li are studied on the basis of the a + A + n + p four-body model. 
We employ the two-body interactions which reproduce the observed properties of any subsystems 
composed of aN, aA and aNN, and aAN. Furthermore, the AiV interaction is adjusted so as to 
reproduce the + -l + splitting of in ^H. The calculated energy splittings of 3/2+-1/2+ and 7/2+-5/2+ 
states in A Li are 0.69 MeV and 0.46 MeV, which are in good agreement with the resent observed 
data. The spin-dependent components of the AA interaction are discussed. 



I. INTRODUCTION 

It is a fundamental problem in hypernuclear physics 
to explore the features of underlying interactions be- 
tween hyperons (Y) and nucleons (N) through analysis of 
many-body phenomena, because YN scattering data in 
free space are quite limited. Then, quantitative analyses 
for light A hypernuclei, where the features of AiV inter- 
actions appear rather straightforwardly in observed level 
structures, are of a special significance. In this connec- 
tion, it is very important that accurate measurements 
for 7-ray spectra have been performed systematically 
0- 0- S Bi bl; which can be used to extract the spin- 
dependent components of AN interactions through the 
detailed analyses of hypernuclear structures. Although 
several shell-model calculations for light A hypernuclei 
have been performed with the restricted configuration of 
(0s) 4 (0p)"0s A USUI!, their structures can be rep- 
resented most excellently by the cluster models. Today, 
it is possible to perform fully microscopic calculations of 
three- and four-cluster systems with sufficient numerical 
accuracy. Such sophisticated calculations make it possi- 
ble to study underlying AiV interactions in comparison 
with the hypernuclear data observed in the 7-ray exper- 
iments 0, 13 ■ Because both short- and long-range cor- 
relations of the A in nuclei are treated very accurately 
in our approach, the characteristics of the free-space AN 
interactions can be studied very precisely. 

The aim of this work is to analyze the ground 1/2+- 
3/2+ and excited 5/2 + -7/2 + doublets in 7 A Li keeping the 
consistency with the 5/2 + -3/2 + doublet in A Be, and to 
determine the AiV spin-spin and spin-orbit interactions 
accurately based on the experimental data for ^Be, A Li 
and A H. The splitting energies for the ground 1/2+-3/2+ 
and excited 5/2 + -7/2 + doublets in A Li 0,0] are related 
intimately to the spin-dependent potentials of the AA 



interaction. Considered naively, the former splitting is 
determined by the spin-spin interaction between the 0s- 
A and the deuteron cluster, while the latter is related to 
both spin-spin and spin-orbit interactions. Thus, it is 
critical whether or not this 5/2 + -7/2 + splitting in A Li 
is reproduced consistently with those of A Li(l/2 + -3/2 + ) 
and A Be(5/2+-3/2+). 

Before starting a realistic calculation with the mi- 
croscopic four-body cluster model, we emphasize that 
the experimental data of the A Li and A Be energy lev- 
els are of a great value for the AN interaction study. 
First the low-lying state energies, A Li(l/2+, 3/2 + ), 
A Li(5/2+,7/2+), A Be(l/2+), and A Be(3/2+, 5/2+), are 
known recently with amazingly high resolution through 
the 7-ray measurements P, 01E HLHl- Secondly, the ba- 
sic structure of these states are well understood on the 
basis of the symmetry consideration without assuming a 
specific form for the underlying AN interactions. 

In order to verify level-energy consistency in the sec- 
ond point mentioned above, let us make a preliminary 
calculation based on the naive SU3 wave functions. In 
other words, we check whether the use of the low-lying 
state energies known for A Li (l/2 + , 3/2 + , 5/2 + ) and A Be 
(1/2+, 3/2+, 5/2+) leads to the right position of the A Li 
(7/2 + ). Based on the nuclear core wave functions 

$ 6 ( 6 Li: l+,3+;T = 0) = | [2] (20) L=0 , 2 (S = 1); J c )(l.l) 
$ 8 ( 8 Be;0+,2+,T = 0) = |[4](40) L=0 , 2 (S = 0); J c ),(1.2) 



the hypernuclear states in A Li and A Be can be assumed to 
have the configurations with an s-state A weak coupling: 
[$jv(L, S; J c , T) <g> A(s!/ 2 )] j+ ■ By using the Hamiltonian 

Ha = Hjsr + £a(si/ 2 ) + ^ Vzva> the hypernuclear level 
energies Ea(Jh) can be expressed straightforwardly as 
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where v a , Vb, v c , Vd, and w e stand for the NA interaction 
matrix elements 

(P3/2S^ /2 |VAr A b3/2^ /2 )2-, {P3/2st/ 2 \VNA\P3/2st /2 )l-, 
(Pl/2S^ /2 \V NA \pi/2sf /2 )l-, (Pl/2Sy 2 \V NA \p 1/2 S^ /2 } - , 

and (p3/ 2 s^ 2 \V N a\p 1 / 2 s^ 2 ) 1 - , respectively. If we input 

the experimental energies for Eq, Eg, Ej, and Eg, a set 
of the interaction matrix elements {v's} are obtained. 
Then, one finds that the use of the solution {v's} leads to 
the theoretical result of 2.53 MeV for E 7 (7/2+) and also 
that this result is quite consistent with the experimental 
value 2.521 ± 0.04 MeV0. Thus these energy levels of 
^Li and ^Be are surely based on the similar maximum 
spatial symmetry, and therefore a detailed realistic cal- 
culation with the microscopic cluster model should have 
an important value. 

In the past, two types of cluster model calculations 
have been performed on the basis of the a + d + A and 
^He+n + p c onfigurations. Using the a + d + A cluster 
model[l(J,[ll|,[l3], it has been discussed that the ground 
1/2+-3/2+ doublet in ^Li is an important candidate to 
extract the AA^ spin-spin interaction. In Ref. 0], we 
proposed the ^He + n + p model, where the full n — 
p correlation is taken into account without the frozen- 
deuteron cluster approximation. The splitting energies 
of the ground and excited doublets in ^Li were calculated 
to be 0.87 MeV and 0.81 MeV, respectively. 

As for the spin-orbit splittings in ^Be and ^ 3 C, the 
three-body (2a + A) model and the four-body (3a + A) 
model have been applied, respectively [l4|. The calcu- 
lated splitting energies of ^Be(5/2 + -3/2 + ) were 80-200 
keV when we adopted the symmetric (SLS) and anti- 
symmetric (ALS) spin-orbit interactions derived from the 
Nijmegen OBE models. These theoretical values are con- 
siderably larger than the experimental value 43 ± 5 keV 
At the same time^H, however, we tried to enlarge 
the ALS potential to be 85 % of the SLS as inspired by the 
a quark-model AN interaction, we predicted the smaller 
splitting of 35-40 keV before the experiment. The similar 
discussion was given also to explain the small spin-orbit 
splitting of the A p-state observed in A 3 C0,0]. 

In this work, we extend these cluster models to the 
four-body treatment of a + n+p+A so as to take account 
of the full correlations among all the constituent parti- 



cles. Such an extended calculation has been tried once in 
Ref. ^5|. Here we focus our attention especially to two 
spin-doublets in ^Li, keeping the consistency with that 
in ^Be. In order to extract the dynamical information on 
the underlying AN spin-spin and spin-orbit interactions, 
two-body interactions among constituent particles (a, n, 
p, A) are chosen so as to reproduce accurately the ob- 
served properties of all the subsystems composed of aN, 
a A and aNN, and a AN. 



II. FOUR-BODY CLUSTER MODEL AND 
METHOD 

In this work, the hypernucleus 7 A Li is considered to be 
composed of a cluster, A particle and two nucleons (N). 
The core nucleus a is considered to be an inert core and 
to have the (0s) 4 configuration, 'J' (a). The Pauli princi- 
ple between the valence nucleon and the core nucleons is 
taken into account by the orthogonality condition model 
(OCM) [3, as the valence nucleon's wavefunction should 
be orthogonal to that the core nucleon. 

Nine set of the Jacobian coordinates of the four-body 
system of ^Li are illustrated in Fig.l in which we further 
take into account the antisymmetrization between the 
two nucleons. The Schrodinger equation are given by 






C=7 



C=8 



C=9 



FIG. 1: Jacobian coordinates for all the rearrangement chan- 
nels (c = 1 ~ 9) of the q + A + n + p four-body model. Two 
nucleons are to be antisymmetrized. 



(H-E)VjmC a U) = , 
H = T + J2 Vab + Vp auli , 

a,b 



(2.1) 
(2.2) 



where T is the kinetic-energy operator and U a b is the in- 
teraction between the constituent particles a and b. The 
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Pauli principle between the a particle and two nucleons 
is taken into account by the Pauli projection operator 
Vpauii which is explained in the next section as well as 
V a b- The total wavefunction is described as a sum of am- 
plitudes of the rearrangement channels Fig. 1 in the LS 
coupling scheme: 



<S>jm ( a Li) 
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Here the operator An stands for antisymmetrization 
between the two nucleons. Xi(-/Vj) and Xi(A) are the 
spin functions of the i-th nucleon and A particle. Follow- 
ing the Gaussian Expansion Method (GEM) [HE El, 
we take the functional form of 4> n im{^)-, i/>nlm(R) and 
&0>c) as 



ip NLM (R.) 



r l e -(r/r 



Y hn {v) , 
R L e -Wn») 2 Y LM (R), 
p * e -(pM 2 Y Xfl (p) , (2.4) 



where the Gaussian range parameters are chosen to lie in 
geometrical progressions: 
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These basis functions have been verified to be suited for 
describing both the short-range correlations and the long- 
range tail behaviors of few-body systems 0, E, H9| . 
The eigenenergy E in Eq.(2.1) and the coefficients C in 
Eq.(2.3) are to be determined by the Rayleigh-Ritz vari- 
ational method. 

The angular-momentum space of the wavefunction 
with l,L,X < 2 was found to be enough for getting satis- 
factorily convergence of the binding energies of the states 
studied below (note that no truncation is taken of the in- 
teractions in the angular-momentum space). As for the 
numbers of the Gaussian basis, n max ,N max and j-Wx, 
4 — 10 are enough. 



III. INTERACTIONS 



The Pauli principle between nucleons belonging to a 
and valence nucleons is taken into account by the orthog- 
onality condition model (OCM) 16]. The OCM projec- 
tion operator Vp au ii is represented by 



Vpauii = lim A \4> 0s (r Na ))(4>(r' Na )\ 

A— »OC 



(3.1) 



which excludes the amplitude of the Pauli forbidden state 
</>0s(r) from the four-body total wavefunction [2l|. The 
Gaussian range parameter b of the single-particle 0s orbit 
in the a particle is taken to be b = 1.358 fm so as to repro- 
duce the size of the a particle. In the actual calculation, 
the strength A for Vp a uli is taken to be 10 5 MeV, which 
is large enough to push the unphysical forbidden states 
into the very high energy region while keeping the phys- 
ical states unchanged. Usefulness of this Pauli operator 
method of OCM has been verified in many cluster-model 
calculations. 



B. NN interaction 

In order to study the fine structure of our a + n+p + A 
system (^Li) , it is necessary to use an NN interaction 
which reproduces accurately the energy spectrum of the 
a + n + p subsystem ( 6 Li). Such an NN interaction 
is given here as follows: We start from the AV'8 [2^ 
potential, Vnn , derived from the AV18 [23| by neglecting 
the (L ■ S) 2 ) term. In our model, this potential gives the 
calculated values of -3.38 MeV and -0.98 MeV for the 
1 + and 3 + states of 6 Li, respectively, being of less binding 
compared to the experimental data. Next, we adjust the 
central and tensor parts of Vjvjv together with the slight 
modification of Vn so that the experimental energies of 
6 Li (l + ,3 + ) and deuteron are reproduced. 



C. qA interaction 

The interaction between the A particle and a cluster 
is derived by folding the AV G-matrix interaction with a 
three-range Gaussian form into the density of the a clus- 
ter in the same manner as our previous work in Ref . |24| . 
In the present work, we employ the G-matrix interaction 
for Nijmegen model F(NF) |2q . the parameters of which 
are also listed in Ref. |2J]. Even if the versions for the 
other Nijmegen models are used, the obtained results are 
almost the same as the present one. This is because our 
AV folding interaction is adjusted so as to reproduce the 
experimental value of £?A( A He). 



A. aN interaction 

For the interaction Vno between a and a valence 
nucleon, we employ the effective potential proposed in 
Ref. which is designed so as to reproduce well the 
low-lying states and low-energy scattering phase shifts of 
the an system. 



D. AV interaction 

For AiV interactions, meson-theoretical models have 
been proposed on the basis of the SU(3) symmetry of 
meson-baryon coupling constants. In principle, these re- 
alistic interactions can be used directly in our four-body 
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model of ^Li. However, the purpose of this work is to ex- 
tract the information on the spin-dependent parts of the 
AN interaction as quantitatively as possible using the 
measured splitting energies of spin-doublet states. We 
employ effective AN single-channel interactions simulat- 
ing the basic features of the Nijmegen meson-theoretical 
models NSC97f [2fi| |. in which some potential parameters 
are adjusted phenomenologically so as to reproduce the 
experimental data. 

Our AN interactions, composed of central, SLS and 
ALS parts are represented as 

^(*-)=EE t, i°' ) exp(-(r/A) 2 ) , (3.2) 

a i—1 



TABLE I: Parameters of the AiV interaction defined in 
Eq.(3.2). Range parameters are in fm and strength are in 
MeV. The numbers in parentheses are improved even state 
strength so as to reproduce observed spin doublet state in 
and odd state strength so as to reproduce observed ground 
doublet state in A Li . 
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Vffi, = V SiS LS+ + U AiS LS_ (3.3) 

with S± =sa±sn- Here, the central potential, V^ N , 
with three-range Gaussian forms are given separately for 
spin-parity states of a= 3 E (triplet even), 1 E (singlet 
even), 3 (triplet odd), 1 (singlet odd). The even- and 

odd-state spin-spin interaction are defined by {V^ N — 

vi^ } )/4 and (V^ 0) - vffl)/*, respectively. 

The potential parameters in the central parts are cho- 
sen so as to simulate AN scattering phase shifts calcu- 
lated by NSC97L The determined parameters are given 
in Table I. It should be noted here that the AN-'EN cou- 
pling interactions are included explicitly in NSC97f, and 
their contributions in many-body systems are different 
from those in free space. This means that our obtained 
phase-shift equivalent potential should be modified ap- 
propriately in applications to hypernuclear system: Wc 

adjust the second-range strengths uij ^ an d v i ^ so that, 
calculated energies of + -l + doublet state in our NNNA 
four-body calculation reproduce the experimental values 
obtained by those of ^H. In Table I, the adjusted val- 
ues of v% E ^ an d v i E ^ are shown in parentheses. On the 
other hand, there was no clear experimental data to de- 
termine quantitatively the odd-state parts, which leads 
to remarkable differences among theoretical interaction 
models. Our present analysis for the splitting energies of 
I Li gives some constraint on the odd-state spin-spin part. 

The second-range values of v% °^ and v% °^ in parenthe- 
ses are determined on the basis of the j^Li data, as shown 
later. 

The SLS and ALS interactions here are chosen so as 
to reproduce the ^Be data. In Ref.Q, the various sets 
were derived from the Nijmegen models. However, the 
5/2 + -3/2 + splitting energies obtained from these sets are 
considerably larger than the experimental value. Now, 
our SLS and ALS interactions are derived as follows: 
Firstly, the SLS part derived from NSC97f with the G- 
matrix procedure is represented in the two-range form 
v sls = j2 2 i=i „(+) e -(xhi) The values of the parame- 
ters are ^ +) = -110.6 MeV and v { 2 +) = -1157. MeV for 
7i = 0.70 fm and 72 = 0.40 fm, respectively as given in 



Ref.Q. Next, assuming V ALS = ~aV SLS , the param- 
eter a is chosen so as to reproduce the measured 5/2+- 
3/2+ splitting energy with the 2a + A cluster model de- 
veloped in Ref.0]. Our obtained value is a — 0.83. Us- 
ing these V SLS and V ALS , we also calculated the energy 
splitting of 1/2--3/2- doublets in ^ 3 C based on 3a + A 
four-body model to be 0.2 MeV which is consisten with 
the observed data within the error This ALS inter- 
action is fairly stronger than that derived from NSC97f. 
As discussed later, however, the similar result can be ob- 
tained by weakening the SLS part without changing the 
ALS part, because only the sum of SLS and ALS is fixed 
by the a Be data. There is a famous quark-model pre- 
diction [23 that the ALS is so strong as to substantially 
cancel the LS one. It should be noted that this prediction 
is not necessarily proved by our present analysis. 



IV. RESULTS 

In Fig. 2, we illustrate our result for the l/2 + -3/2 + and 
5/2+-7/2+ doublet states of ^Li • The energies of the 1 + - 
3 + doublet state of 6 Li nucleus calculated in the frame- 
work of the a+n+p three-body model are —3.7 MeV and 
— 1.6 MeV, being measured from the a+n+p three-body 
threshold. As shown in the left side of the figure, the cal- 
culated splitting energies for both doublets are about 1 
MeV very similar to that for + -l + doublet state of ^H 
(^He), when only the even-state central interaction is 
used. Namely, the even-state spin-spin interaction turns 
out to contribute similarly to the + -l + splitting energy 
of (^He) and the 1/2+-3/2+ and 5/2+-7/2+ ones of 
ILL 

Next, let us switch on the odd-state central interac- 
tion. When only the even-state interaction is used, the 
obtained value of the ground-state energy is —9.79 MeV. 
When we use the 1 and 3 interactions derived from 
NSC97f, the ground 1/2+ state is obtained at -9.23 
MeV. This energy changes only slightly (—0.06 MeV) 
with the inclusion of SLS and ALS, because the spin- 
orbit interactions have essentially no effect on the 1/2+ 
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FIG. 2: Calculated energy levels of ^Li on the basis of Q + A + 
n+p model. The energies are measured from the a + A + n+p 
threshold. The observed energy splittings of 3/2 + -l/2 + and 
7/2 + -5/2+ are 0.69 MeV and 0.47 MeV, respectively. 



state due to its L = structure. The final value —9.29 
MeV means that the experimental A binding energy (5.58 
MeV) is reproduced well, because the calculated energy 
of the a+p+n subunit is —3.7 MeV in our model. Found 
in the figure, this fact is owing to the peculiar role of our 
odd-state interaction. Therefore, the repulsive contribu- 



tion from the spin-independent part (SV^^ + V^')/A 
turns out to be decisive to reproduce the experimental 
value. The repulsive nature of this part is an important 
property of NSC97f. On the other hand, the correspond- 
ing part in the recent model ESC04 [28| is attractive in 
contrast to NSC97f. The important role of the repulsive 
odd-state interaction in our analysis does not necessarily 
mean that the odd-state part in NSC97f is more realistic 
than the one in ESC04, as discussed later. 

As for the l/2 + -3/2 + splitting, the addition of the 
NSC97f odd-state central interaction leads to 0.97 MeV. 
This splitting is too large in view of the experimen- 
tal value (0.69 MeV), because the contribution of the 
SLS/ALS interactions to this l/2 + -3/2 + splitting is quite 
small. We add here the attractive (repulsive) correction 
on the 3 ( 1 0) state interaction, which works efficiently 
in high-spin (low-spin) partners of doublets, by making 
the spin-spin interaction more attractive: We introduce 
the attractive spin-spin interaction in second range (0.8 

fm), Av ss = -100.0 MeV. The modified values of v% 0) 
( 3 o) 

and v 2 are shown in parentheses of Tabled, which leads 
to the calculated values of 0.71 MeV and 0.54 MeV for 
the lower and higher doublets, respectively (cf. Fig. 2), 
as seen in the figure. 

Now, we come to the important stage of looking at 
the roles of the SLS and ALS interactions for splitting 



energies. It should be noted here that these interactions 
work differently for the two doublet states of ^Li: The 
contributions to the ground-state 1/2+-3/2+ doublet are 
very small, where the pn pair part outside the a core 
is dominated by the L — component spatially. On 
the other hand, in the case of the excited 5/2 + -7/2 + 
doublet composed of the L = 2 pn pair, the SLS and ALS 
interactions play important roles: As seen in Fig. 2, the 
SLS works attractively (slightly repulsively) for the 7/2+ 
(5/2+) state, because the 7/2+ state is dominated by the 
spin-triplet configuration of the L — 2 pn pair and the 
s-state A. On the other hand, the ALS works efficiently 
in the 5/2 + state which has both configurations of spin- 
triplet and spin-singlet. The ALS which acts between 
5 = and S = 1 AN two-body states has essentially no 
effect on the 7/2+ state. 

Thus, owing to the combined effects of the SLS and 
ALS, our final result reproduces nicely the observed en- 
ergies of the spin-doublet states in ^Li. Recently, Millncr 
calculate 3/2+-1/2+ and 7/2 + -5/2 + spin-doublets states 
using a shell model with the phenomenological interac- 
tion matrix element 0, His calculation is also in 
good agreement with the recent data. It should be noted 
here that the strength of the ALS part is not necessar- 
ily determined by our present analysis. The above result 
is obtained by making the ALS part stronger than that 
given by NSC97f so as to reproduce the ^Be data. The 
very similar result, however, can be obtained by weaken- 
ing the SLS part without changing the ALS part. 

Before going to summary, we would like to comment on 
the role of the A7V-S7V coupling. Our basic assumption 
in this work is that the AN-HN coupling interaction can 
be renormalized into the AN-AN interaction effectively. 
In this spirit, the even-state parts of our AN-AN inter- 
action were adjusted so as to reproduce the + and 1 + of 
^H. It is reasonable, however, to consider that the AN- 
HN coupling works more repulsively in ^Li. It is likely 
that the role of the odd-state repulsion in our treatment 
is a substitute for this effect. This is the reason why the 
attractive odd-state interaction in ESC04 models can- 
not be denied. As shown in Fig. 2, the energy of 5/2+ 
state is located above by about 0.2 MeV in comparison 
with the observed energy of 5/2+ state. This problem 
may be solved by taking into account the repulsive effect 
of the AN-T,N coupling instead of the odd-state repu- 
sion, because the SLS/ALS interaction works more effi- 
ciently under the attractive odd-state interactions. Some 
authors |30l l3l| pointed out the extra contribution to 
the ^H(0 + -1 + ) splitting energy from the three-body cor- 
related AN-HN mixing. The present authors also ob- 
tained the value of about 0.3 MeV for the three-body 
contribution of AN — E7V coupling in the + -l + split- 
ting energy in |29j . In the shell model calculation 
|32j |. Millener calculated the spin-doublets states in ^Li 
including AN — E7V coupling and he concluded that this 
contribution was small in these splitting energies. On the 
other hand, Fetisov pointed out that the large effect of 
AN — S N coupling was found in both of jji and 7 A Li [3j| . 
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It is an open problem to study A7V-E7V coupling effects 
consistently for and ^Li. 

In summary, we discussed the two spin-doublets of 
3/2+-1/2+ and 7/2+-5/2+ in ^Li based ona + A + n+p 
four-body model. Here, it is important that all the two- 
body interactions are chosen so as to reproduce both the 
binding energy of any subsystem composed of two- and 
three-constituent particles. Our AN interactions, simu- 
lating AY scattering phase shifts calculated by NSC97f, 
are adjusted so as to reproduce the observed data of spin- 
doublet states. It is found that the even-state AN inter- 
action leads to the similar values of the splitting energies 
of the 0+-1+ doublet in QHe) and the 1/2+-3/2+ and 
5/2+-7/2+ doublets in ^Li. Then, the odd-state inter- 
actions play important roles to reproduce the difference 
between the two doublet states in ^Li. With use of the 
SLS and ALS interactions adjusted so as to reproduce 
the 5/2+-3/2+ splitting in ^Bc, the two doublet states 



in ^Li can be reproduced exactly by tuning the odd-state 
spin-spin interaction. 

The basic assumption in our present approach is that 
the AN-Y,N coupling interactions are renormalized rea- 
sonably into our AN interactions. The validity of this as- 
sumption will be investigated in our future studies. The 
coupled four-body calculation of a + A + N + N and 
a + E + Y + AMsin progress. 
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